Hadron-quark phase transition in asymmetric matter with boson condensation 
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In the present work we study the hadron-quark phase transition with boson condensation in 
asymmetric matter by investigating the binodal surface and extending it to finite temperature in 
order to mimic the QCD phase diagram. We consider a system with two conserved charges (isospin 
and baryon densities) using the Gibbs' criteria for phase equilibrium. In order to obtain these 
conditions we use two different models for the two possible phases, namely the non-linear Walecka 
model (NLWM) for the hadron matter (also including hyperons) and the MIT bag model for the 
quark phase. It is shown that the phase transition is very sensitive to the density dependence of 
the equation of state and the symmetry energy. For isospin asymmetry of 0.2 and a mixed phase 
with a fraction of 20% of quarks, a transition density in the interval 2po < pt < 4po was obtained 
for temperatures 30 < T < 65 MeV. 

PACS numbers: 21.65.-f„25.75.Nq,05.70.Fh,12.38.Mh 



I. INTRODUCTION 

Since some decades ago just after the discovery of 
asymptotic freedom of QCD [1] the possibility of the ex- 
istence of a new state of matter in high energy physics 
is under consideration, namely, a color deconfined phase 
of quarks and gluons, the so-called quark-gluon plasma 
(QGP) [2]. The main goal of the heavy-ion collision ex- 
periments at ultra-rclativistic energies is to create, under 
controlled conditions, and understand the properties of 
this new state of matter. This opens a new field of study 
in strong interaction physics. 

Many experiments have been proposed and acceler- 
ators built in the search for the QGP at different en- 
ergies at SIS/GSI, AGS/BNL, SPS/CERN, RHIC/BNL 
and LHC / CERN to look for some signs and signatures of 
the production of the QGP that subsequently hadronizes 
[3]. The study of particle production in ion collisions 
contribute to the understanding of the conditions under 
which the quark-gluon plasma may be produced and also 
to determine the equations of state (EoS) of strongly in- 
teracting matter. 

In hydrodynamical models the system which arises 
from a high-energy collision (fireball) reaches an approx- 
imately local thermal equilibrium (thermalization) and 
expands evolving collectively up to the point when the 
mean free path of the created and interacting particles 
becomes large enough for the particles to escape from the 
fluid, i.e., the interactions among the particles of the sys- 
tem cease because the system has reached the freeze out 
point. The approximately local thermalization is con- 
sidered due to detailed computations of the expansion 
stage that takes much longer than the typical scattering 
times [4]. Although, non-equilibrium processes are also 
important for the dynamics, equilibrium processes are a 
quite good approximation to be used in theoretical mod- 
els, and are a reasonable first approximation at freeze- 
out. Some authors consider both the temperature at 



which the inelastic collisions cease (chemical freeze-out) 
and the elastic collisions cease (kinetic freeze-out) [5]. 

In the end of the eighties of the last century the ex- 
istence of the (chiral) critical end point (CEP) in the 
QCD phase diagram was suggested [6, 7] and since then 
its properties have been extensively studied [8]. Although, 
most lattice QCD calculations indicate the existence of 
the CEP for fi B > 160 MeV [9-11], its exact location 
is not well known since it depends, for example, on the 
mass of the strange quark. The CEP separates second- 
order transition at high temperatures (or even a smooth 
cross-over) from the first-order transition at high chemi- 
cal potentials in the QCD phase diagram. Studying this 
intermediate region is a hard task since perturbation the- 
ory cannot be applied to QCD at this regime and ad- 
ditionally at finite chemical potential the usual lattice 
approach fails. Moreover, new techniques have been pro- 
posed to study lattice QCD at finite T and /1 [12]. On 
the other hand the lattice QCD simulations of different 
groups disagree with each other on the location of the 
critical end point. 

Subsequently in the late nineties a hypothesis arose 
[13] that the onset of the deconfinement phase transi- 
tion was located between the top AGS and SPS ener- 
gies. The CERN energy scan program of the NA49 ex- 
periment at SPS has given signs of a phase change at 
E\ a h ~ 30 A-GcV particularly from the horn-like peak in 
the K+/tt+ ratio [14]. 

Furthermore, the hadronic freeze-out estimated for 
different colliding energies [15] shows a maximum at 
■y/sjvjv =4 + 4 GeV, which can be reached for a fixed- 
target bombarding energy of 20 — 30 A-GcV at the bary- 
onic chemical potential region /i# = 400 — 500 MeV. In 
addition, hydrodynamical calculations [16, 17] of phase 
trajectories during collisions, in the QCD phase diagram, 
indicate that for E uh ~ 30 A-GeV {y/s^ ~ 8 GeV) the 
trajectory goes near the CEP. 

Since then, the interest on the intermediate energies 
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TABLE I: Ion beam top energies in some collision experiments. 



(not ultra-relativistic) in collision experiments is increas- 
ing as well as the theoretical study of the phase transition 
at that regime. 

For this purpose the new facilities, namely, the 
Nuclotron-based Ion Collider f Acility (NIC A) at JINR/ 
Dubna [18] and the Facility for Antiproton and Ion 
Research (FAIR) at GSI/Darmstadt [19], give the op- 
portunity to explore an interesting region of the phase 
diagram, in search for the QGP and where the criti- 
cal end point (CEP) is expected to exist, complement- 
ing other heavy-ion collision experiments (NA61-SHINE, 
low-energy RHIC) compatible with the energy range 
E\ a b = 2 — 40 A- GeV. It is also important to consider the 
possibility that new features arise at still lower temper- 
atures and higher densities as the color superconducting 
quark matter phases like the 2SC phase and the recently 
conjectured quarkyonic phase [20]. A summary of ion 
beam top energies used in some collision experiments is 
shown in table I. 

It is also possible to study the phase transition from 
hadronic matter to a quark phase within the effective 
models that describe two separated phases, and also the 
structure of the mixed phase can be obtained through 
the Gibbs' conditions [21]. Some features of this phase 
transition can be obtained by means of the binodal sur- 
face, which is a phase coexistence curve in the parameter 
space. 

The different ion beams used in collision experiments 
present different numbers of neutrons (N) and protons 
(Z). It is also interesting to study the isospin effects on 
the transition to a mixed phase of hadrons and quarks. 
We can define the asymmetry parameter (isospin ratio) of 
a nucleus (or the hadron phase) as a = (N — Z)j (N + Z). 
such that, a runs from (symmetric matter) to 1 (pure 
neutron matter) . From table II one sees some ions used in 
nucleus-nucleus collisions and the respective asymmetry 
parameter of each system. Systems with isospin ratios 
< a < 0.23 are up to now experimentally accessible 
in ion collisions and the case a = 1.0 corresponds to 
neutron matter which is relevant in some astrophysical 
applications. 

We study the phase transition from hadrons to a 
quark-gluon plasma in asymmetric matter using a two- 
phase model, analyzing the features that depend on the 
isospin and may be relevant in a phenomenological de- 



scription of heavy-ion collisions [22-25] . 

It is interesting to investigate asymmetric systems 
since in the liquid-gas phase transition of nuclear matter 
the asymmetric case shows different properties from the 
symmetric one [22, 26, 27]. It is shown that the transition 
of an asymmetric system is of second order (continuous) 
rather than the first order (discontinuous) transition in 
symmetric systems [22-25]. 

Hence, an interesting task is to investigate the isospin 
effect on the hadron-quark phase transition at lower tem- 
peratures and densities higher than the saturation den- 
sity of the normal nuclear matter, that can be probed 
in heavy-ion collisions at intermediate energies. In ad- 
diction, the presence of bosons can modify the isospin of 
the hadron phase. Also, at low temperatures these fea- 
tures depend strongly on the nuclear symmetry energy. 
On the other hand, at higher temperatures the inclusion 
of bosons shows an interesting feature due to the onset of 
a boson condensate in asymmetric systems if we consider 
an approximately local thermal equilibrium. 
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TABLE II: Some ions used in collision experiments and the 
respective asymmetry parameter (a) of the system. 



This approach is useful for providing a qualitative ori- 
entation on the features that arise when a phase transi- 
tion from hadrons to quarks takes place and two con- 
served charges are considered, i.e. at finite baryon den- 
sity and isospin. 

As already mentioned, the problem we investigate 
in the present paper has already been studied in pre- 
vious works [23-25] within different perspectives, based 
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on different paramctrizations and containing different in- 
gredients. In [23] the hadronic phase is given by one 
paramctrization of the non-linear Walccka model, the 
quark phase is calculated with the MIT bag model for 
one specific value of the bag constant and pions are in- 
cluded. In [24, 25] a mixed phase of hadrons and quarks 
is particularly emphasized and the influence of the sym- 
metry energy on the phase transition investigated. In 
[25] neither bosons nor gluons were considered and the 
quark phase was described within the MIT bag model, 
and in [24] two quark different models have been used: 
the MIT bag model (with and without gluons) and the 
color dielectric quark model. In these works five differ- 
ent parametrizations of the non-linear Walecka model are 
considered and one of them includes the 5 mesons. All 
these models have a quite high value of the symmetry 
energy slope, namely 85 < L < 103 McV, and therefore 
they have a quite hard symmetry energy at intermediate 
densites, the densities of interest for the present work. 

In the present work we consider seven different 
parametrizations of the non-linear Walecka model, which 
span a large variety of EOS: they include hard and soft 
EOS with hard and soft symmetry energies at intermedi- 
ate densities. This will allow us the see the effect of both 
the isoscalar and the isovector interaction on the phase 
transition. For the quark model we have considered the 
MIT bag model with various values of the bag constant 
and with gluons. The bag constant was chosen in ac- 
cordance with heavy-ion collision data. In the hadronic 
phase we have studied the effect of including two kinds of 
bosons, pions and kaons. This more complete picture al- 
lows us to discuss various aspects of the phase transition 
at finite temperature which were not discussed before. 

The remainder of this article is organized as follows: 
In section II we present the formalism used in this work. 
In section III the mixed phase features are presented and 
in section IV we show the numerical results and discus- 
sion. Finally, in section V we summarize the results and 
give a brief concluding discussion. 



Quark matter has been extensively described by the 
MIT bag model [28]. In its simplest form, the quarks 
are considered to be free inside a Bag and the thermody- 
namic properties are derived from the Fermi gas model 



in two limits: T = 0, m„ ^ and T =/= 0, 



0. The 



energy density, the pressure and the quark q density are 
respectively given by: 



£ a = 3 x 2 



q—u^d 



d 3 p 



V 2 +m 2 q (f q+ + f q -) + B. 



P q = 



dp 



(2.1) 



(/,+ + /,-) -fl, (2.2) 



3x2 



d 3 p 
(2^p 



(/<+-/*-) 



fq± = J 



1 



+ e Ke„TI* q )/T] 



(2.3) 



(2.4) 



where 3 stands for the number of colors, 2 for the spin 
degeneracy, m q for the quark masses, B represents the 
bag pressure and f q ± the distribution functions for the 

quarks and anti-quarks, e q = \/p q +m^, ±/x g being the 



chemical potential for quarks and anti-quarks of type q, 



(2.5) 



The quark density is 



n u + n d , 



(2.6) 



II. THE FORMALISM 



and the "quark baryon density" is given by: 



In the present section we present the equations of 
state (EoS) for the hadron phase and for the quark phase 
used in this work and their respective definitions. Bosons 
are included using a meson-exchange type Lagrangian 
that couples the bosons to meson fields and the possibil- 
ity of a boson condensate is also presented. 

A. Quark phase: quarks u, d (+ gluons) 



rid 



(2.7) 



The thermodynamic potential per unit volume of the 
MIT bag model (two-flavor case) and the corresponding 
equations of state (EoS) [23, 29] for massless quarks and 
a Bosc gas of gluons of degeneracy 7 g = 2 x 8 with the 
lowest-order gluon interaction (a s ) is 
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from where we can obtain the pressure P Q 



-Qqgp/V the energy density and the quark number density: 



P, 
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' 45" T I 1 "17 



15a, 
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2 T ^ + — 



-B 



(2.9) 



£qgp — 3-Pq G p + 4i? 



(2.10) 



The strong coupling a s is taken as a constant in the 
present work (a s = 0.349) and Nf stands for the number 
of flavors (Nf = 2, quarks u and d). 

B. Hadron phase: nucleons (+ hyperons) 



The equations of state of asymmetric matter within 
the framework of the relativistic non-linear Walecka 
model (NLWM) [30] are presented next. In this model 
the nucleons are coupled to neutral scalar er, isoscalar- 
vector w M and isovector-vector p^ meson fields. The La- 
grangian density reads 



O-B 1 2 2 , 1 i 3 , 1 \ 4 1 2 2 1 ^ 4 



E/ (27r )3 { 



In 



l + e-ZJ^-'i) 



In 



+ e-n'!^]}-Kg l ,M&, (2 . 12) 



where /3 = 1/T, E) = (p 2 + M* 2 ) 1 / 2 and the effective 
chemical potential of baryon j is given by 



+ -d.ad^a - -mlcr 2 - -ka 3 - -Xa* 
2 2 3! 4! 



1 



1 



(2-11) 



K7 = Mi - fc w o - t 3J g p P03 • (2.13) 
The EoS for the baryons can then be calculated 



= —Y 

3tt 2 ^ 



p A dp 



3 - x /p* + m 

,2 



= {fFj+ + fFj-) + -^^0 
#2 ^ 



£ 4 , m p 2 m o- 2 ^3 

24^° + ~T Po3 ~ ~ a ° ~ 6 a ° 



24 ff0 



where m* = m 3 — g a j a is the baryon effective mass, 

O/ii/ = dp^v - , = - - g p (p,j, x 
<7ij are the coupling constants of mesons i = a, w, p with 
baryon j, and m, is the mass of meson i. The couplings 
fc (fc = 2 Mjv 5 3 an d ^ — 6 9^ c ) are the weights of 
the non-linear scalar terms and r is the isospin operator. 
This Lagrangian includes an isoscalar-isovector mixing 
term A v (g 2 p p . p^)(g 2 uj^uj^) as presented in [31] which 
plays an important role in high densities. It can also 
be extended to include all the hyperons from the baryon 
octet. 

Within the relativistic mean field (RMF) framework 
the thermodynamic potential per unit volume corre- 
sponding to the Lagrangian density (2.11) is 



+ A v g 2 p gl Wq P03 j 



(2.14) 



£b = ^E / P 2 dpy/p 2 + m?(fF j + + If,-) 



£ 4^ m 'P 2 
g w + — P03 "+ 

OA 2 2 2 2 
3Av 9 p 3a; w P03 



k 3 , A 4 

6 24 



(2tt) 3 



(/fj+ - fp 3 -)d 3 p 



(2-15) 



(2.16) 
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and /fj± is the Fermi distribution for the baryon (+) 
and anti-baryon (-) j: 



fFj± = 



(2.17) 



C. Hadron phase: bosons (pions + kaons) 

It is possible to include the boson fields using terms 
from the chiral perturbation theory [33]. In the present 
work we prefer to use a meson-exchange type Lagrangian 
that couples the bosons to meson fields. For simplicity 
we apply the same approach to the kaons and pions. The 
Lagrangian density in the minimal coupling scheme [34- 
39] is given by: 



C b = D* Df<S>- mT $* $ 



where the covariant derivative is: 



d u + iX, 



X, 



9ub + 9pb U ■ Pfj, 



(2.18) 



(2.19) 



(2.20) 



and the boson effective mass, ml = m b — g ab a. The 
boson field can then represent either the kaons or pions 



J 



(particles and anti-particles): 

$ = (K + ,K a ) , $* = (K-,K°), (2.21) 



or 



$=(7T-,7r°) , $* EE (7T + ,7r°) . (2.22) 

The isospin third-component to the bosons is given by 



T"3tt = 




1~3K = 



(2.23) 



In order to obtain the boson thermodynamic potential it 
is also possible to perform a similar calculation as carried 
out in reference [39] with the respective modifications in 
the covariant derivative (2.19). 

As the neutral pion is its own anti-particle we need 
to set Xfj_ = [g w „o = g pn o = 0, and if required 
to* o = m,o , then g a ^o = 0) to achieve the correct ther- 
modynamical features of these uncharged particles. In 
this case, the Lagrangian (2.18) takes up its simplest form 
and $* = <f>. In particular, ir° results completely decou- 
pled from the other particles and its population is that 
of a boson gas at temperature T and /i^o = 0. 

In the Appendix we show the calculation of the 
bosonic EoS: 



= C 2 [^ b -X )^m^]-T j 7^{ln 



1 _ e -/3(^ + -M) 



In 



1 _ e -P(" +m) 



(2.24) 



£ b = ( 2 [ml 2 + nl - X 2 )] + 



(2.25) 



n b 



2C 2 (fib -Xo) + j ^3 { f B+ - f B -} , (2.26) 
where the Bose distribution for particles (Jb+) and anti-particles (Jb~) appears naturally in the EoS and read: 



e/ 9(^±TMb) _ 1 e /3[(e*±X )T«,] _ l e^Tn) _ 1 ' 



with el = \fp 2 + ml 2 , and hence we define the boson effective chemical potential as 

Vb = ^b — X . 



(2.27) 



(2.28) 



From equation (2.26) one notes two contributions in the boson density and we can define them as the "condensate" 
and "thermal" ones 



n b = n c b (()+n b r (T) 



(2.29) 
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and the entropy density is given by s b = /3(P b + £b — Pb n b) ■ The order parameter £ can be obtained through the 
minimization of the thermodynamic potential. 



D. Hadron phase equations 

The thermodynamic potential of the hadron phase (HP) 
including both the baryons and the bosons, is given by 

n HP = n B + n b (2.30) 

where flp is given by (2.12) and fl b by (A19). By min- 
imizing the thermodynamic potential £Ihp with respect 
to the meson fields er, uj and p, and also with respect to 
the order parameter £, within the mean-field approxima- 
tion (<r (a) = (T ; ^ = Sfto ujq ; p^ ->• (p M ) = 
<5m0 ^* 3 Po = ^m0^ i3 P03)j w e obtain the equations for the 
hadron phase: 



depending on the signal of p b (either positive or nega- 
tive). Thus 

p b = ml + X Q or p b = -{ml - X ) , (2.37) 
and 

p b - X = ml or p b - X = -ml , (2.38) 

such that the condition for the onset of the condensate 
state is 

v b — > ml or Vb — > —ml . (2.39) 



^2gab{n c b + n s b ) 



ml w 



tot 

6 



6 Hfe 



OA 2 2 2 

- 2A v 5 p 3 W P03 w o , 

m p z' 03 = T3 j n j + X! T3fc n;> 

OA 2 2 2 

- 2A V g p 5w w po3 , 



and 



C [pb-^b(o)] [pb + u b 



where 



d 3 p 



(2^) 3 E, 



HfF+ + fF-), 



(2.31) 
(2.32) 

(2.33) 



is the baryon scalar density of particle j, and the respec- 
tive baryon density 



d 3 p(f F+ -f F -) 



(2.34) 



1 ~ (27T) 3 

The "boson scalar density" for the boson b is given by 

^=jfCfi» +/.-), (2.35) 

and the boson density is given by (2.26). 

From the last equation of (2.32) we obtain the condi- 
tions for the possibility of a boson condensate (( = 0, no 
condensate) , resulting 

p b = uj^(p = 0) or p b = -uJ b (p = 0) , (2.36) 



According to (2.24) and (2.39) the condensate (zero mo- 
mentum state) does not contribute to the pressure of the 
system as expected. When the condensate is not present, 

C = o. 



III. THE MIXED PHASE 

In the following, three situations for the phase coexis- 
tence are discussed in detail: A) hadron matter consti- 
tuted of nucleons and quark matter constituted of quarks 
u and d, B) hadron matter constituted of nucleons and 
pions and quark matter constituted of quarks u and d 
and C) hadron matter constituted of nucleons, hyperons, 
pions and kaons with zero net strangeness and quark mat- 
ter constituted of quarks u and d. 



A. A. Nucleons and quarks 

According to the Gibbs' conditions [21] for the phase co- 
existence the chemical potentials, temperatures and pres- 
sures have to be identical in both phases (H = hadron 
phase; Q = quark phase): 



u Q 



J<H _ rpQ 

P H (pX,pX,T) = pQ(pQ,plT) 



(3.1) 



The conservation of the isospin (n 3 ) and baryon densi- 
ties (tib) are also required, so that in terms of these two 
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charges [22] and including the mixed phase we can write 
P H (n%,n*,T) = pQ(n%nlT), 

t,»(n*,n»,T) = »<i(n%nlT), 
^(n*,n*,T) = $(n%nlT), (3.2) 

n B = (1 - %)nf + x^b i 
n 3 = (1 - x)nf + , 

where for the hadron phase 



where n T = n^. + and we assume <7 p jv = g P = g-n- 
In order to obtain the simplest thermodynamic features 
for the pions we set g^ — and for simplicity g alx = 
so that in this case m* = m„. The in-medium (s-wave) 
Bose effective pion energy is 



Wtt- (p = 0) = m ff - g p P03 



(3.9) 



and the pion chemical potentials and the effective tt 
chemical potential 



(3.10) 



and for the quark phase 



q _ n u - n d 



"3 



(3.3) 



(3.4) 



|U w o = 



(3.11) 



As jU n > /i p then /i^- > and according to (2.39) the 
onset of the pion (tt~ ) condensation takes place when 



v-- -> m. 



and the EoS for the hadronic phase 



Ph =Pe 



P. 



£h — £i 



(3.12) 



(3.13) 



and x represents the fraction of quarks in the mixed 
phase. The asymmetry parameter, a (isospin ratio) of 
the nuclei was defined as 



where P^ and £ v are given by (2.24) and (2.25) and we 
have also included the neutral pions as a free Bose gas. 



N — Z n n — n p 
N + Z ~ n~B 



(3.5) 



and the asymmetry parameter of the hadron and quark 
phases can be defined by: 



■2 ^ ; a« = -2^, (3.6) 



,9 



hence 



H n n n p q nd — n u 
a = ■ ; a w = 3 ■ , (3.7) 

n n +n p n d + n u 



such that < a H < 1 (just nucleons case) and the quark 
one < a Q < 3. 



B. Nucleons, pions and quarks 

When bosons are present the isospin density of the 
hadron phase is modified according to (2.31) and a H can 
be greater than 1. The isospin density of the hadron 
phase with ir~ becomes 



,H 



(3.8) 



C. The baryon octet, pions, kaons and quarks 

At this stage we have included in the hadron phase 
all baryons of the baryon octet and in order to keep 
the strangeness conservation as J^. Si = in both phases 
we also have included the K + meson in the hadron 
phase. For the Gibbs' conditions (3.1) we need to add: 
|Ug = jUy , so that we can write the chemical potential 
as fii = Bi/j,B + hi A*3 + Sif-is, where Bi^I^i and Si are 
the baryonic, isospin and strangeness quantum numbers 
of particle i. 

The equations for the baryons and bosons are already 
presented in this work. For the kaons we set g p x = g P , 
guK = and also g a K = such that m* K = ttlk as for 
the pions. We are aware that this choice is very naive. It 
was done in order to explore the isospin degree of free- 
dom. In a future work a more realistic parametrization 
of the kaon-meson coupling will be used which will al- 
low us to discuss the strangeness degree of freedom more 
completely. 



IV. RESULTS 

First of all it is important to present some features 
of the MIT bag model. Figure 1 shows a qualitative 
overview of the MIT bag model in a simple case when 
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FIG. 1: (a) A simple qualitative overview on the T x \ib 
curve in the MIT bag model (dashed curves) for the case 
Hu = , compared with the freeze-out curve (continuous 
line) from figure (b). The bag constant values from I to IV are: 
B 1/4 = 145, 160, 190 and 210 (MeV). (b) A parametrization 
of the freeze-out curve deduced from particle multiplicities in 
heavy- ion collisions, Cleymans et al. (2006) [40]. 



Mu = P"d i comparing with data analysis of some colli- 
sion experiments [40]. This figure also indicates the best 
values of the bag constant B to be used in some energy 
ranges that we describe in the present work. In high tem- 
peratures and low baryon chemical potential (low den- 
sity) we use the values B 1 / 4 = 190 and 210 MeV. On 
the other hand, our analysis at intermediate energies is 
performed with B 1 / 4 = 160 MeV. 

For the hadronic phase we use the parameter sets pre- 
sented in table III, where we give the symmetric nuclear 
matter properties at saturation density as well as the 
parameters of the models. In Figs. 2 (a) and (b) the 
pressure of symmetric nuclear matter and the symmetry 
energy, respectively, are plotted for a large range of densi- 
ties. In Fig. 2 (a) we also include the experimental con- 



straints obtained from collective flow data in heavy-ion 
collisions [46] . We have considered a wide range of mod- 
els frequently used to study stellar matter or finite nuclei. 
Even though some of them do not satisfy the constraints 
determined in [46] , as a whole these sets of models allows 
us to understand the influence of a hard/soft equation 
of state (EOS) and a hard/soft symmetry energy of the 
hadron matter-quark matter phase transition. We have 
considered: NL3 [44] with a quite large symmetry energy 
and incompressibility at saturation and which was fitted 
in order to reproduce the ground state properties of both 
stable and unstable nuclei, TM1 [41] which also repro- 
duces the ground state properties of both stable and un- 
stable nuclei and provides an equation of state of nuclear 
matter similar to the one obtained in the RBHF (Rela- 
tivistic Brueckner Hartree-Fock) theory, softer than NL3 
at high densities, TMlojp [42], the TM1 parametriza- 
tion with a mixed isoscalar-isovector coupling which we 
fix in order to obtain a softer density dependence of the 
symmetry energy [31], FSU [32] which was accurately 
calibrated to simultaneously describe the GMR in 90 Zr 
and 208 Pb, and the IVGDR in 208 Pb and still reproduce 
ground-state observables of stable and unstable nuclei; 
GM1 and GM3 [45] generally used to describe stellar mat- 
ter, with a symmetry energy not so hard as the one of 
NL3 and TM1, and NLp [43], which has been used to dis- 
cuss the hadron matter-quark matter transition in [25], 
and which presents an EOS at high densities between 
GM1 and GM3. 

Let us first describe some hadron-quark matter sys- 
tems at zero and finite temperature, including the decon- 
fined phase transition, through isothermal processes. We 
first discuss the effects of pions and gluons on the phase 
transition. For this discussion we take NL3 to describe 
the hadronic matter, however the main conclusions do 
not depend on the nuclear model considered. 

In Figs. 3 (a) and (b) we show slices of the binodal 
surface indicating the two-dimensional phase-coexistence 
boundary in {P, T, a} space, at T = 30 MeV. For each 
temperature, the binodal section is divided into two 
branches. One branch describes the system in the hadron 
phase, while the other branch describes the quark-gluon 
phase. In both figures the role of the pions and gluons is 
presented. Fig. 3 (a) shows a system with no gluons for 
two cases with and without the pions. The blue curves 
represent a system with no pions. The gluons have a 
very strong effect on the critical point, corresponding to 
a maximum in the pressure, when both phases coexist for 
a = 0. The presence of gluons increases the critical pres- 
sure almost by 100%. The role of the pions is better seen 
in Fig. 3 (b) (red lines) . The asymmetry parameter of the 
hadron phase increases due to the presence of the pions 
which increase the isospin interaction. In equilibrium, 
the pressures in both phases must be equal according 
to the Gibbs' conditions. When pions are present these 
conditions still hold. We observe a slight increase of the 
pressure of the quark-gluon phase for a > 3, following 
the hadronic pressure increase. 
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FSU [32] 


TM1 [41] 


TMlup [42] 


NLp [43] 


NL3 [44] 


GM1 [45] 


GM3 [45] 


n (fm 3 ) 


0.148 


0.145 


0.145 


0.160 


0.148 


0.153 


0.153 


K (MeV) 


230 


281 


281 


240 


271.76 


300 


240 


m* /m 


0.62 


0.643 


0.643 


0.75 


0.60 


0.70 


0.78 


m (MeV) 


939 


938 


938 


939 


939 


938 


938 


-B/A (MeV) 


16.3 


16.3 


16.3 


16.0 


16.299 


16.3 


16.3 


£ sym (MeV) 


32.6 


36.9 


31.9 


30.5 


37.4 


32.5 


32.5 


L (MeV) 


61 


110 


55 


85 


118 


94 


90 


m CT (MeV) 


491.5 


511.198 


511.198 


512 


508.194 


512 


512 


m u (MeV) 


782.5 


783 


783 


783 


783 


783 


783 


m p (MeV) 


763 


770 


770 


763 


763 


770 


770 




10.592 


10.029 


10.029 


8.340 


10.217 


8.910 


8.175 




14.302 


12.614 


12.614 


9.238 


12.868 


10.610 


8.712 




11.767 


9.264 


11.147 


7.538 


8.948 


8.196 


8.259 


b 


0.000756 


-0.001506 


-0.001506 


0.006935 


0.002052 


0.002947 


0.008659 


c 


0.003960 


0.000061 


0.000061 


-0.004800 


-0.002651 


-0.001070 


-0.002421 




0.06 


0.0169 


0.0169 














Av 


0.03 





0.03 















TABLE III: Parameter sets used in this work and corresponding saturation properties. 



At finite temperature pions are present as a Bose gas 
and their presence as a condensate state at low enough 
temperatures is also possible. The presence of a pion gas 
and a pion condensate changes the pressure at low densi- 
ties according to Fig. 4 by increasing the absolute value 
of the p meson field (i.e., the isospin interaction) since 
the condensate itself does not contribute to the pressure 
of the system as a boson gas. The lowest pressures of the 
binodal occur for the largest values of the asymmetry 
parameter a (1 for the hadronic phase without pions). 

When gluons are included in the quark phase the den- 
sities reached by the system at the binodal surface in- 
crease slightly in both phases such that the onset of the 
pion condensation takes place at a slightly higher den- 



sity: 2.87no instead of 2.18no at T = 30 MeV according 
to Fig. 5, when an isothermal process is analyzed. There- 
fore, the presence of gluons shifts the phase transition to 
a quark-gluon plasma to larger densities. 

The onset of the pion condensation according to the 
cq. (2.39) [or similarly eq. (2.36)] is clearly seen in Fig. 6 
where we plot the pion mass m„, pion chemical poten- 
tial fj, w , pion effective chemical potential v„ and pion fre- 
quency at p = 0, u>£(p = 0). The pion condensation 
occurs for the lower densities when the conditions (2.39) 
or (2.36) are satisfied. 

In Fig. 7 we show the binodal slices at different temper- 
atures and for the bag constant 190 MeV. The enclosed 
area becomes smaller with increasing temperature and 
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FIG. 3: Binodal sections at T = 30 MeV and the effect of pions and gluons. The NL3 parameter set is used. The blue lines 
indicate a system with no pions. The red ones indicate the presence of pions. (a) Without gluons. (b) With gluons. 
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FIG. 4: Binodal section: pressure as a function of the bary- 
onic density. The presence of the pions increases the isospin 
interaction such that a pressure increase is observed. 
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FIG. 5: Population of particles at T = 30 MeV with and 
without the gluons. The filled circles mark the onset of the 
pion condensate at 2.18no in the first case and at 2.87no in 
the second one. 



the pressure at a = decreases when the temperature in- 
creases. The two branches merge into a single line when 
the system reaches the critical temperature at zero chem- 
ical potential and density. The critical temperature (T c ) 
of the phase transition is ~ 150 MeV for the bag con- 
stant B 1 / 4 = 190 MeV. For larger values of B we obtain 
a larger pressures at the same temperature and the other 
way round for smaller values. The results shown in the 
figure are consistent with the ones found in ref. [23] al- 
though here the NL3 parameter set has been used. The 
calculation includes both pions and gluons. 

Next we discuss the inclusion of strangeness. The 
population of particles at T = 50 MeV can be seen in 
Fig. 8 for two cases: (a) a simple system of protons, p, 
neutrons, n, and pions, 7r~, in the hadron phase, and 
(b) including the hyperons of the baryon octet and K + 



mesons in the hadron phase. In both cases the total 
strangeness of the system is zero, therefore, we just have 
quarks u and d in the quark phase. Fig. 8 (a) shows an 
increase of pions at low baryon densities, which plays an 
important role in the isospin density of the system. Most 
of the pions below 2.6nrj are in a zero momentum state 
(i.e., a pion condensate). The same pattern can be seen in 
Fig. 8 (b) on pions and nucleons, indicating that strange 
particles are not important in these conditions at that 
temperature but they do appear at higher densities. We 
do not see kaon condensation, just a pion condensate as 
in the first case. It is important to analyse how sensitive 
are the above results to the choice of the kaon-meson 
interaction. Work in this direction will be done in the 
near future. 
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FIG. 6: The onset of the pion condensation for the case with 
gluons in Fig. 5. The y-axis we plot the pion mass, chemical 
potential, effective chemical potential and frequency at p = 0. 
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FIG. 7: Slices of the binodal surface using the NL3 parameter 
set and B 1/4 = 190 MeV. The critical temperature is T c ~ 
150 MeV. The calculation includes both pions and gluons. 



We now discuss the effect of the density dependence of 
the EOS on the binodal surfaces. In Fig. 9 (a) one sees a 
comparison of the hadron phase-quark phase binodal sec- 
tions among the different parameter sets listed in Table 
III for the zero temperature case and B 1 / 4 = 160 MeV. 
Qualitatively all the curves behave in the same way. The 
difference lies in the pressures and the densities reached 
by the different systems which is explicitly shown in 
Fig. 9 (b). We conclude that the different behaviors seen 
for the binodal sections are due to the EoS at large den- 
sities, see Fig. 2 (a) where the pressure is plotted as a 
function of density for cold symmetric nuclear matter. 
At finite temperature a similar trend is obtained except 
that the maximum densities reached are smaller. 

The effect of the symmetry energy on the binodal is 
better discussed analyzing Fig. 10 where the binodal for 
TM1 and TMlwp is plotted without pions. These two 



models have the same isoscalar behavior and just differ 
in the isovector channel, TMlwp having a softer symme- 
try energy. We conclude that a softer symmetry energy 
favors a phase transition at larger asymmetries. For the 
same reason the models with a softer symmetry energy 
have their binodals for larger asymmetries in Fig. 9 (a). 

Fig. 2 (a) and (b) are used in the following to dis- 
cuss the differences between the models. We see that the 
hadron density at the binodal surface is very sensitive 
to the softness/hardness of the EoS at intermcdiatc/high 
densities. In particular, the largest pressures are attained 
by the softest EOS. It is interesting to analyse the behav- 
ior of TM1: it behaves at low densities as a hard EOS 
like NL3 and at high densities as a soft one, giving the 
largest pressure at the critical point. It is the relative 
change of hard/soft character of the EOS that explains 
the crossing between the different models in Fig. 9 (a). 
We have not included a curve for FSU because due to 
its softness no phase transition was obtained at reason- 
able densities. The behavior at large densities can be 
adjusted by changing the value of the parameter x which 
multiplies the forth power of the w-meson term in the 
Lagrangian density. A larger value gives a softer EOS at 
large densities. We have reduced the value of \ and for 
X = 0.03 we could get convergence at reasonable densi- 
ties. This coincides with the large density behavior of 
the new parametrization proposed in [49], that corrects 
the behavior of FSU at large densities which predicted 
too small maximum star masses and too large star radii. 

The density dependence of the energy density docs 
not affect the binodal surface of symmetric nuclear mat- 
ter but it certainly has an effect if we consider asym- 
metric matter. We investigate the phase transition at 
intermediate energies using a convenient choice of differ- 
ent parametrizations of the NLWM in order to explore 
different compressibilities at large densities as well as an 
asysoft and asyhard EoS. We take into account the pa- 
rameter sets: NL3, hard EoS and symmetry energy; NLp, 
intermediate behavior both in the isoscalar and isovector 
channel; TM1, soft EoS at high densities and hard sym- 
metry energy and TMlwp, with a soft symmetry energy. 

In order to discuss the effect of isospin asymmetry 
on the binodal sections, we allow the temperature to 
change with fixed asymmetry parameter, and compare 
the predictions of the different models. Figs. 11 (a) 
and (b) show, for NL3 and NLp, the binodal sections in 
{ng, T, a} space and the projection of several branches at 
different a onto the (ng,T) plane (HP = hadron phase; 
QP = quark phase). In other words Figs. 11 (a) and 
(b) show the QCD phase diagram with different asym- 
metries, a = 0, 0.2, 0.4, 0.6, 0.8, and 1.0, from the right 
(I) to the left (II) in the two phases. From now on, in or- 
der not to reach too high densities in the hadron phase we 
exclude the gluons from the system. This does not affect 
the comparison between models and may give rise to a 
maximum 20% underestimation of the transition density. 

In Fig. 11 (a) we have considered B 1 / 4 = 190 MeV 
together with the models NL3 and NLp. The properties 
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FIG. 8: Population of particles for the case in Fig. 7 (a) at T = 50 MeV. (a) Just p, n and n in the hadronic phase, (b) The 
baryon octet, ir~ and K + in the hadronic phase. 




of the EoS are clearly reflected on these results: for NL3 
the transition occurs for smaller densities due to its very 
large compressibility at large densities. It is also this high 
value of the compressibility that dilutes in part the effect 
of the asymmetry parameter. NLp has a much softer EOS 
and symmetry energy and, therefore, the curves obtained 
for a fixed asymmetry span a larger range of densities. In 
summary, the hadron-quark phase transition is favored 
when the asymmetry of the system is increased. 

We are interested in discussing the phase transition 
at intermediate temperatures and high densities and, for 
this reason, we consider smaller bag pressures according 
to Fig. 1. We set B 1 ^ 4 = 160 MeV in order to reach 
a specific range of temperature and densities, which is 
presented in Figs. 12 (a) and (c) and also in Fig. 13 (a). 
The asymmetries experimentally available up to now ac- 



cording to table II are in the range — 0.23. 

Since NL3 is too hard and does not satisfy most of the 
constraints imposed by experimental and observational 
measurements [47] we consider in the following TM1 with 
and without a cop non linear term which allows us to 
discuss a asy-soft and a asy-hard EoS. We also take into 
account the NLp parametrization in order to compared 
with the results already obtained in [24]. 

In Fig. 12 we compare TM1 and TMl-wp. This allows 
us to discuss the effect of density dependence of the sym- 
metry energy on the phase transition since the isoscalar 
channel is kept fixed. The main effect of a softer sym- 
metry energy is to shift the binodal sections for larger 
values of the asymmetry parameter to larger densities. 
A harder symmetry energy allows the occurrence of the 
hadron-quark phase transition at smaller densities, and 
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therefore, easier to reach with heavy ion collisions at in- 
termediate energies. Similar conclusion were drawn in 
[25] where the effect of the (5-meson on the phase tran- 
sition was discussed: the <5-meson gives rise to a harder 
symmetry energy at large densities favoring the hadron- 
quark phase transition. 

In Fig. 13 (a) we show for the same bag constant the 
binodal sections obtained with NLp. It is seen that due to 
a softer EoS at intermediate densities the binodal sections 
occur at larger densities when compared with TM1. The 
effect of the bag constant is clear if we compare Fig. 11 
(b) with B 1 / 4 = 190 MeV with Fig. 13 (a). A larger 
B shifts the phase transition to much larger densities, 
showing that in order to obtain a good estimation it is 
essential to choose an adequate value of B. 

We show in Figs. 12 (b) and (d) and Fig. 13 (b) a part 
of Figs. 12 (a) and (c) and also Fig. 13 (a), corresponding 
to a = 0.2. We also include curves corresponding to the 
mixed phase with the quark concentrations \ = 0.2 and 
0.5 which correspond to 20% and 50% of quarks in the 
mixed phase. One sees the indication of an interesting re- 
gion where the phase transition probably occurs and can 
be probed by intermediate energy heavy-ion collisions. 
This region is located in the range ns = 2 — 4 no and 
T = 50 — 65 MeV and can reached by the new planned 
facilities (NICA) at JINR/ Dubna [18] and (FAIR) at 
GSI/Darmstadt [19] that will start operations in the next 
few years. 

The density behavior at intermediate/high densities 
defines the transition region. For instance, models TM1 
and TMlwp would favor the detection of a quark phase 
more than NLp. 



V. SUMMARY 

We have presented a study of the deconfinement phase 
transition from hadronic matter to a quark-gluon plasma, 
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FIG. 11: (a) Binodal section in (ns, T, a) space and the 
projection of several branches for different asymmetry pa- 
rameter a onto the (ns, T) plane. The asymmetries are: 
a = 0, 0.2, 0.4, 0.6, 0.8, 1.0, from the right (I) to the left 
(II) in the two phases, with no gluons. (b) Same as figure (a) 
for the NLp parameter set. In both cases the critical temper- 
ature where p,B — is T c ~ 155 MeV. 



which could be formed in heavy-ion collisions. Calcu- 
lations at finite temperature with a simple two-phase 
model and the inclusion of pion and kaon condensa- 
tion were done in order to describe this type of system. 
We have studied the effect of the density dependence of 
the EoS on the phase transition choosing a convenient 
set of parametrizations of the NLWM. We have consid- 
ered both hard and soft EoS at intermediate densities 
as well as models with asyhard and asysoft symmetry 
energies. We have also considered the effect of gluons 
on the quark phase. For the quark phase we have used 
the MIT bag model and chose the bag constant accord- 
ing to a parametrization of the freeze-out curve deduced 
from particle multiplicities in heavy-ion collisions [15]: 
for deconfinement phase transition at T ~ 50 — 60 and 
p ~ 2 — 6po: the bag constant B 1 / 4 ~ 160 MeV was used. 
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FIG. 12: Same as in Fig. 11 for B 1/4 = 160 MeV. The critical temperatures are T c ~ 130 MeV and the labels (I) and (II) also 
represent the asymmetry in the same way as in Fig. 11. (a) The TM1 parameter set is used, (b) Part of Fig. 12 (a) for a — 0.2 
with the mixed phase for different quark concentrations (x = 0.2, 0.5). (c) The TM1 parameter set and the mixing term A v 
has been used in the present case, (d) Part of Fig. 12 (c) with the mixed phase. 



An important result is the difference between the 
phase diagram for a symmetric system and that for asym- 
metric matter as observed in liquid-gas phase transition. 
Usually, the onset of the phase transition takes place at 
lower baryon densities and temperatures in more asym- 
metric systems. This can be probed by means of neutron- 
rich nuclei in heavy-ion collisions. Moreover, the density 
at which the phase transition occurs is sensitive to the 
density dependence of the EoS at intermediate densities. 
A hard EoS gives rise to a transition at lower densities. 
The density dependence of the symmetry energy also af- 
fects the transition when asymmetric matter is consid- 
ered. The phase transition is favored for asymmetric nu- 
clear matter, and even more for an asyhard symmetry 
energy. 

Both thermal pions and pion condensation have been 
included in the calculation. They mainly play a role at 



low densities, large isospin asymmetries and large tem- 
peratures. We have considered that the pions couple to 
the nucleons through the p-meson [23] . Using an equiva- 
lent parametrization for the kaon-meson coupling, which 
maybe too naive since it only takes into account the 
isospin interaction, we have verified that the effect of in- 
cluding strangeness in the hadron phase was negligible for 
a system with an overall strangeness equal to zero. This 
can be generalized to finite strangeness when it becomes 
possible to prepare heavy-ion collisions with hypcrnuclci. 
It remains to be investigated how sensitive are the results 
to the pion and kaon interaction. 

The results obtained for the phase transition are very 
sensitive to the EOS. Both the isoscalar and isovector 
interactions have an effect on the transition density. Ac- 
cording to the effective models used in this work there 
exists a region in the parameter space where the phase 
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transition probably occurs and can be probed by heavy- 
ion collisions at intermediate energies. This region is lo- 
cated in the range Ub = 2 — 4 no and T — 30 — 65 MeV 
and can be reached by the new planned facilities (NICA) 
at JINR/ Dubna [18] and (FAIR) at GSI/Darmstadt [19] 
that will start operations in the next few years. We have 
obtained a larger T interval, extending to lower temper- 
atures, for the same densities obtained in [25] due to the 
properties of the models used. A clear sign of a phase 
transition could be used to constrain both the EOS and 
symmetry energy at intermediate densities. 

We have verified that models with a soft EOS and 
soft symmetry energy such as FSU do not predict a 
hadron-quark phase transition at densities that could be 
attainned in the laboratory. 

A more complete system with all baryons of the bary- 
onic octet and strange mesons, as well as interacting pi- 
ons and kaons, and using interactions constrained by ex- 
perimental measurements, is under investigation in order 
get more systematic results. 



it is possible to obtain the respective thermodynamic po- 
tential and the EoS of the boson fields. It is convenient 
to transform $ into real and imaginary parts using two 
real fields, (f>i(pc,t) and 02 (x, t) such that 



$ = — (fa+ifa) 



^ = -^{cj> x -icj> 2 ). (A2) 



The conjugate momenta are 



d(d <f>i) 



(A3) 
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Appendix A 
The boson thermodynamic potential 

Using the Lagrangian density in the minimal coupling 
scheme [34-39] 



and the corresponding Hamiltonian density of the boson 
field, Tib = irido4>i + 7r 2 <9o02 — C-b such that the four- 
current and its zero component are 

j„=i [**(D M *)-(DJ**)$] , (A4) 

.70 = <j>2Kl ~ 017T2 ■ (A5) 

For the neutral pions we just have ir = = dn4> an d 

$ = i sucn that = $ and j M = 0. Now we can 
write the Hamiltonian density 
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n b = + \y 2 + ni(Xo <h) - <k) + \{^<t>i? + ^(V0 2 ) 2 + {di ( j)2)x l 0i 



. t 2 



(di fa)x t 4>2 + M fa) 2 + M fa) 2 + -%-(<t>l + 4>l) , 



(A6) 



where i = 1,2,3, and the partition function in the grand canonical ensemble as a functional integral is given by 



Z b = J [d-Ki][d-K 2 ] J [dfa][dfa] exp | J dr J 



d 3 o 



*7Ti-5 1" l7r 2-^ fib + ^b(fa^\ - 0itt 2j 



<9t 



<9t 



(A7) 



where is the boson chemical potential associated with the conserved charge Q = J d 3 x jo(x). Here "periodic" means 
that the integration over the field is constrained in the imaginary time variable t = it so that 0fc(x, 0) = 0fc(x, 
and where /? = 1/T. The neutral pion Hamiltonian is H n o = ^ir 2 + i(V0) 2 + |m^ 2 which has the form of that of 
a neutral scalar field, so that it can be used within the relativistic mean field approach, as it is known that the pion 
pseudoscalar interaction term vanishes in the mean field level. After some algebra the integration over momenta can 
be done and the result is 



[dfa][dfa] exp{ I dr I d 3 x 



d0 
dr 



- - i (n b - X ) fa 



d<i>2 



i (lib - X ) 0i 



(A8) 



\{Vfaf - ^(V0 2 ) 2 + (9,0i) A ?; 2 - (a i 2 )X i 0i - \(X t fa) 2 - \(X t fa) 2 - \m* b 2 (fai + 2 ) 



where N is a normalization factor. In the mean field approach (Xi) = 0. Integrating (A8) by parts, and taking into 
account the periodicity of 0i and 2 , the result is 



Z b =N 2 



[#i][# 2 ] expli / dr I d 3 x{ 



dr 2 



V 2 - ml 2 + (jib - X f 



— + V 2 -m* 2 + ( Mfc -A ) 2 



2 + 2i(p b - X ) 



dr 



dr 



(A9) 



The fields can be expanded in Fourier series as 



1/2 

MX,T) =V2C COS(0) + ( £ ) ^5> i;(P ' X+ "" T) 01,n(p) 



(A10) 



where the Matsubara frequency is u n = 2irnT, due to the constraint of periodicity of the fields, such that 0fc(x, /3) = 
0/c(x, 0) for all x. The normalization factors of (A10) can be chosen so that each Fourier amplitude is dimensionlcss. 
The infrared character of the field is carried out by ( and 6, so that, 0i,o(p = 0) = 2 ,o(p = 0) = which allows some 
particles to reside in the n = 0, p = state, i.e., a possibility of a condensation of the bosons into the zero-momentum 
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state ("s-wave" condensation). Using (A10) in (A9), and noting that (/>_„(— p) = c/>*(p) because </>i(x,r) and </>2(x, r) 
are real fields, we have 



Z h = N z 



where 



n p 



[ / #l >Tl (p)#2,n(p) 



5 = /3UC 2 [(//6-Xo) 2 -m 6 * 2 ] -2 EE [^i,-n(-P) > <h,-n(~P) 

n p 



01,n(p) 



(AH) 



(A12) 



and 



D = f3 2 



U)l + p 2 + m* b 2 - Of, - X ) 2 -20*6 - A )w„ 



(A13) 



2{fi b - X Q )uj n wl + p 2 + m* b 2 - (M6 - A ) 2 

As the thermodynamic potential is given by O = — (1//3) ln(Z), we can perform the integrals in (All) and write 



\n(Z b ) = PVC 2 [(p b - X ) 2 - m* b 2 ] + In [(det D)"* 



(A14) 



The multiplication of by any constant is irrelevant since it does not change the thermodynamics of the system. 
The second term of (A14) is given by 



iln[detD] = -iln|nil^[( 

I n p 



( w 2 + p 2 + m* b 2 - - X ) 2 Y + 4( Mb - X ) 2 lo 2 



- \ In ( J] /? 2 K + (- + - ^) 2 ] } - \ ^ { II H + (' 



(A15) 



,. - -. zi 

so that (A14) can be written as 

ln(Z b ) = PVC 2 [0*6 " A ) 2 - m£ 2 ] - 1 ]T In {/3 2 [ W 2 + ( u + Mb ) 2 ] } - \ E ^ {/? 2 K + + M & ) 2 ] } , (A16) 

n,p n,p 

and in the continuum limit, neglecting the zero-point energy contribution, due to the mean field approach, the result 
is 



In(Z 6 ) = PVC 2 [0x 6 - A ) 2 - m£ 2 ] - U | ^ { In 



1-e" 



In 



1-e" 



(u> +fi) 



where 



ui±lp) = Jp 2 + mf ± A 



(A17) 
(A18) 



is the effective Bose energy, such that the thermodynamic potential for the bosons is given by 



fib _ ln(Z b ) 

v ~ W~ 



e K 2 -(M6-A„) 2 ] + r/^{ 



(27T) = 



In 



1 - e 



-/3(w+- M ) 



In 



1 - e 



} ■ (A19) 



[1] D. J. Gross and F. Wilczek, Phys. Rev. Lett. 30, 1343 [2] J. C. Collins and M. J. Perry, Phys. Rev. Lett. 34, 1353 
(1973). 



18 



(1975); B. Freedman and L. McLerran, Phys. Rev. D 16, 
1130 (1977); Phys. Rev. D 16, 1147 (1977); Phys. Rev. D [22 
16, 1169 (1977); V. Baluni, Phys. Lett. B72, 381 (1978); [23 
V. Baluni, Phys. Rev. D 17, 2092 (1978); S. Chin, Phys. [24 
Lett. B78, 522 (1978); J. Kapusta, Nucl. Phys B148, 
461 (1978); E. V. Shuryak, Phys. Rep. 61, 71 (1980); [25 
[3] Saeed-uddin, Eur. Phys. J. C 6,355 (1999); P. Braun- 

Munzinger, I. Heppe and J. Stachel, Phys. Lett. B 465, [26 
15 (1999); P. Braun-Munzinger, D. Magestro, K. Redlich 
and J. Stachel, Phys. Lett. B 518, 41 (2001); D. Zschi- [27 
esche, S. Schramm, J. Schaffner-Bielich, H. Stocker and [28 
W. Greiner, Phys. Lett. B 547, 7 (2002); F. Becattini, 
M. Gazdzicki, A. Keranen, J. Manninen and R. Stock - [29 
Phys. Rev. C 69, 024905 (2004); and many others. 
L. McLenan, Rev. Mod. Phys. 58, 1021 (1986). [30 
B. Mohanty, Nucl. Phys. A830, 899c (2009); B.I. Abelev 
et ai, STAR Collaboration, Phys. Rev. C 79, 34909 
(2009) and references therein. [31 
M. Asakawa, K. Yazaki, Nucl. Phys. A504, 668 (1989). 
A. Barducci, et al, Phys. Lett. B231, 463 (1989); A. [32 
Barducci, R. Casalbuoni, S. DeCurtis, R. Gatto and G. 
Pettini, Phys. Rev. D 41, 1610 (1990); A. Barducci, R. 
Casalbuoni, G. Pettini and R. Gatto, Phys. Rev. D 49, [33 
426 (1994). 

M. A. Stephanov, Prog. Theor. Phys. Suppl. 153, 139 [34 
(2004) [hep-ph/0402115vl]; Int. J. Mod. Phys. A 20, 
4385 (2005); R. Casalbuoni, hep-ph/0610179, and ref- 
erences therein. [35 
Z. Fodor and S.D. Katz, J. High Energy Phys. 04, 050 
(2004) [arXiv:hep-lat/0402006vl]. [36 
R. V. Gavai, S. Gupta, Phys. Rev. D 78, 114503 (2008); 
E. Shintani, S. Aoki and Y. Kuramashi, Phys. Rev. D [37 
78, 14503 (2008); R. V. Gavai, S. Gupta, Phys. Rev. D 
71, 114014 (2005). [38 
Anyi Li et al, arXiv:0908.1155. 

Z. Fodor and S.D. Katz, Phys. Lett. B534, 87 (2002) [39 
[hep-lat/0104001]; Z. Fodor and S.D. Katz, J. High En- 
ergy Phys. 03, 014 (2002) [hep-lat/0106002]. [40 
M. Gazdzicki and M. I. Gorenstein, Acta Phys. Polon. B 
30, 2705 (1999) and references therein; M. I. Gorenstein, 
Phys. Atom. Nucl. 71, 1594 (2008). [41 
M. Gazdzicki, J. Phys. G 30, S701 (2004); S. Afanasiev 
et al, Phys. Rev. C 66, 054902 (2002). [42 
J. Randrup and J. Cleymans, Phys. Rev. C 74, 047901 

(2006) ; arXiv:0905.2824. [43 
Y. B. Ivanov, V. N. Russkikh, and V. D. Toneev, Phys. 
Rev. C 73, 044904 (2006); nucl-th/0503088. [44 
A. S. Khvorostukhin, V. V. Skokov, V. D. Toneev, and K. 
Redlich, Eur. Phys. J. C 48, 531 (2006); nucl-th/0605069. [45 
A. N. Sissakian, A. S. Sorin, V. D. Toneev, nucl- 
th/0608032. See also: http://nica.jinr.ru . [46 
S. Chattopadhyay, J. Phys. G 35, 104027 (2008); P. 
Senger et ai, J. Phys. G 36, 064037 (2009); Johann [47 
M. Heuser for CBM Collaboration, arXiv:0907.2136. See [48 
also: http://www.gsi.de/fair . 

L. McLerran and R. D. Pisarski, Nucl. Phys. A796, 83 [49] 

(2007) ; L. McLerran and R. D. Pisarski arXiv:0706.2191; 
L. McLerran arXiv:0907.4489. 

L. D. Landau and L. Lifshitz, Statistical Physics, Perga- 



mon Press, Oxford, 1969. 

H. Miiller and B. D. Serot, Phys. Rev. C 52, 2072 (1995). 

H. Miiller, Nucl. Phys. A 618, 349 (1997). 

M. Di Toro, T. Gaitanos, V. Greco and A. Lavagno, Nucl. 
Phys. A775, 102 (2006). 

M. Di Toro, B. Liu, V. Greco, V. Baran, M. Colonna, S. 

Plumari, arXiv:0909.3247v2 [nucl-th]. 

M. Barranco and J. R. Buchler, Phys. Rev. C 22, 1729 

(1980). 

N. K. Glendenning, Phys. Rev. D 46, 1274 (1992). 

A. Chodos, R.L. Jaffe, K. Johnson, C.B. Thorne and V.F. 
Weisskopf, Phys. Rev. D 9 (1974) 3471. 

N. K. Glendenning, Compact Stars, Springer- Verlag, 
New- York, 2000. 

B. D. Serot and J.D. Walecka, Adv. Nucl. Phys. 16 (1986) 

I. ; J. Boguta and A. R. Bodmer, Nucl. Phys. A 292, 413 
(1977). 

C. J. Horowitz and J. Piekarewicz, Phys. Rev. Lett. 86, 
5647 (2001). 

B. G. Todd-Rutel and J. Piekarewicz, Phys. Rev. Lett. 
95, 122501 (2005); F. J. Fattoyev and J. Piekarewicz, 
Phys. Rev. C 82, 025805 (2010). 

A. E. Nelson and D. B. Kaplan, Nucl. Phys. A479, 285c 
(1988). 

R. Knorren, M. Prakash, and P. J. Ellis, Phys. Rev. C 
52, 3470 (1995); V. Thorsson, M. Prakash and J. M. 
Lattimer, Nucl. Phys. A572, 693 (1994). 
J. Schaffner and I. N. Mishustin, Phys. Rev. C 53, 1416 

(1996) . 

M. Prakash, I. Bombaci, M. Prakash, P. J. Ellis, J.M. 
Lattimer and R. Knorren, Phys. Rep. 280, 1 (1997). 
N.K. Glendenning and J. Schaffner-Bielich, Phys. Rev. C 
60, 025803 (1999). 

J. A. Pons, S. Reddy, P.J. Ellis, M. Prakash and J.M. 
Lattimer, Phys. Rev. C 62, 035803 (2000). 
Joseph I. Kapusta, Finite- Temperature Field Theory, 
Cambridge University Press, 1989. 

J. Cleymans, H. Oeschler, K. Redlich and S. Wheaton, 
Phys. Rev. C 73, 034905 (2006); J. Phys. G 32, S165 
(2006). 

K. Sumiyoshi, H. Kuwabara and H. Toki, Nuclear Physics 
A581, 725 (1995). 

H. Pais, A. Santos, L. Brito and C. Providencia, Phys. 
Rev. C 82, 025801 (2010). 

B. Liu, V. Greco, V. Baran, M. Colonna, and M. Di Toro, 
Phys. Rev. C 65, 045201 (2002). 

G. A. Lalazissis, J. Konig, P. Ring, Phys. Rev. C 55, 540 

(1997) . 

N. K. Glendenning and S. A. Moszkowski, Phys. Rev. 
Lett. 67, 2414 (1991). 

P. Danielewicz, R. Lacey and W. G. Lynch, Science 298, 

1592 (2002); arXiv:nucl-th/0208016v2. 

T. Klahn, et al, Phys. Rev. C 74, 035802 (2006). 

C. J. Horowitz and J. Piekarewicz, Phys. Rev. Lett. 86, 
5647 (2001). 

F. J. Fattoyev, C. J. Horowitz, J. Piekarewicz, G. Shen, 
arXiv:1008.3030vl [nucl-th]. 



